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ABSTRACT 

The binary system PSR B1259— 63/LS 2883 consists of a 47.8 ms radio pulsar that orbits the com- 
panion Be star with a period of 3.4 years in a highly eccentric orbit. The system has been well sampled 
in radio, X-ray, and TeV 7-ray bands, and shows orbital phase-dependent variability in all observed fre- 
quencies. Here we report on the discovery of >100 MeV 7-rays from PSR B1259— 63/LS 2883 through 
the 2010 periastron passage. Using data collected with the Large Area Telescope aboard Fermi from 
33 days before periastron to 75 days after periastron, PSR B1259— 63/LS 2883 was detected at a sig- 
nificance of 13.6 standard deviations. The 7-ray light curve was highly variable over the above period, 
with changing photon index that correlates with the 7-ray fiux. In particular, two major flares that 
occur after the periastron passage were observed. The onset of 7-ray emission occurs close to, but not 
at the same orbital phases as, the two disk passages that occur ~1 month before and '^l month after 
the periastron passage. The fact that the GeV orbital light curve is different from that of the X-ray 
and TeV light curves strongly suggests that GeV 7-ray emission originates from a different component. 
We speculate that the observed GeV ffares may be resulting from Doppler boosting effects. 
Subject headings: gamma rays: stars — gamma rays: observations — Pulsars: individual 
(PSR B1259-63) — X-rays: binaries 



1. INTRODUCTION 

The binary system PSR B1259-63/LS 2883 com- 
prises a young radio pulsar with a period 47.8 ms 
and a Be star LS 2883. With an eccentric (e~0.87) 
orbit, the pulsa r approaches the periastron every 3.4 
years () Johnston et al.lH994l ). The high spin-down power 
of PS R B1259-63 (~8xl035 erg s"^: IManchester et al.l 
II995I ) suffices to generate a relativistic pulsar wind 
(PW). The system is highly variable over an orbital pe- 
riod in radio (e.g., Johnston et al. 2005), X-ray s (e.g ., 
IHiravama eFaLl |l999; Chc rnvakova et al. 2006, ' 20091 ) . 
and TeV -Y-ravs ([Aharonian et al.l |2005., 2009). The 
broadband electromagnetic spectrum is believed to result 
from the interaction of the PW and the stellar wind of 
LS 2883, the latter being composed of a polar wind and 
a dense equatorial circumstellar disk. The st ellar disk 
is inclined with respect to the orbital plane (Wex et al.l 
119981: 1 Johnston et al.|[r999D such that the pulsar passes 
through the disk shortly before and shortly after the pe- 
riastron passage. 

The X-ray flux changes significantly with orbital phase. 
The 1-10 keV fiux increases from ^^10""'^^ erg cm~^ s~^ 
at apastron to more than '-^lO"^^ erg cm~^ s~^ 
shortl y before and shortly a fter the periastron pas- 
sage ( Chernvakova et all 120091 and references therein). 
The X-ray photon index (Tx) also shows variability re- 
lated to the orbital period and a spectral state with 
the hardest spectra (Fx ~ 1-2) occurred around the 
same time as the observed ra pid growth of the X-ray 
flux (jChernvakova et all I2006D . The unpulsed radio 
flux increases through the past four periastron passages 
with respect to other orbital phases, although the light 
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curves show slightly different behaviors during each pas- 
sage (Johnston et al. 2005, and references therein). 

The first detection of PSR B1259-63/LS 2883 in 7- 
rays was made by the H.E.S. S. Cherenkov arr ay through 
the 2004 periastron pa ssage (lAharoniaii et al.| [2005). and 
subsequently in 2007 (|Aharonian et al.ll2009fl . The 2004 
and 2007 data show that the TeV emission peaks ~10 
days before and ^20 days after the periastron passages 
and the possible dips seen in the TeV light curves before 
and after the periastron se em to coincide with the stellar 
disk passage (see Fig. 1 in lKerschhaggl|[20TTl as well as 
Fig. [II). 

An upper limit of ^10~^° erg cm~^ s~ ^ was derived us- 
ing EGRET data (iHartman et al.|[T999[) for the average 
7-ray emission. A claim made by the AGILE collabora- 
tion on the detection of PSR B1259-63/L S 2883 above 
100 M eV in a 2-day period in August 2010 (jTavani et al.l 
I2OIOD has not been confirmed by data collected using 
the mo re sensitive Fermi /L AT detector during the same 
period (|Abdo et al.ll2010cD . 

During a regular monitoring of 7-ray emission from 
PSR B1259— 63 shortly before the periastron passage in 
mid-December, we found the first evidence (^4ct) for 7- 
ray emission from PSR B125 9-63 /LS 2883 during a 3- 
day time interval (Tam et al.l 120101). T he discovery was 
later confirmed by Ab do et al.l ()2010dl ) who used data 
collected over 30 days from 2010 November 18. 

In this Letter, detailed Fermi analysis results of 
PSR B1259-63/LS 2883 from 2008 through early 2011 
are presented. 

2. FERMI/LAT OBSERVATIONS AND RESULTS 

The Large Area Telescope (LAT) aboard the Fermi 
Gamma-ray Space Telescope can detect 7-rays with en- 
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ergies between ~20 MeV and >300 GeV (jAtwood et all 
l2009f ). The 7-ray data used in this work were ob- 
tained between 2008 August 4 and 2011 February 28. 
These data are available at the Fermi Science Sup- 
port CenteiQ. Due to the discovery of 7-rays from 
PSR B1259— 63/LS 2883, a modified sky survey was com- 
menced from 2010 December 27 for 10 days. In this 
mode the southern hemisphere receives 30% extra expo- 
sure. We used the Fermi Science Tools v9rl8p6 pack- 
age to reduce and analyze the data in the vicinity of 
PSR B1259-63/LS 2883. Only events that are classified 
as the "diffuse" class or the "data-clean" class were used. 
To reduce the contamination from Earth albedo 7-rays, 
we excluded events with zenith angles greater than 105°. 
The instrument response functions "P6_V3_DIFFUSE" 
were used. 

We carried out unbinned maximum-likelihood analyzes 
(gtlike) of the circular region with a 15° radius cen- 
tered on the 7-ray position of PSR B1259— 63 (see be- 
low) . We subtracted the background contribution by in- 
cluding the Galactic diffuse model (glLiem_v02.fit) and 
the isotropic background (isotropic_iem_v02.txt), as well 
as all sources in the first Fermi/LAT catalog (IFGL; 
lAbdo et al.ll2010"bh within the circular region of 25° ra- 
dius centered on the 7-ray position of PSR B1259— 63. 
We assumed single power laws for all IFGL sources con- 
sidered, except for 7-ray pulsars of which the spectra 
follow power laws with exponential cut-off ()Abdo et al.l 
l2010aD . The spectral parameter values for sources within 
10° from PSR B1259— 63 as well as the normalization pa- 
rameters of the diffuse components were set free. 

No significant emission was found before 2010 Novem- 
ber 11. We derived upper limits of 7-ray flux from 
PSR B1259-63/LS 2883 using data obtained between 
2008 August 4 and 2010 November 10 (see Fig. H]). 
However, the source became active roughly when 
PSR B1259-63 en tered the stellar disk (|Tam et al.ll2010l: 
lAbdo et alllMoa . 

We first analyzed 0.2-100 GeV data from 2010 Novem- 
ber 11 to 2011 February 28. Using a power-law 
description of PSR B1259— 63, the ma ximized test- 
statistic (TS) value (jMattox et al.l I1996D obtained for 
the pulsar position is 184, corresponding to a detec- 
tion significance of 13.6tT. This allows us to classify 
PSR B1259-63/LS 2883 as a new GeV 7-ray source. 
The best-fit position of the 7-ray emission is estimated 
by gtfindsrc to be at right ascension (J2000) = 195?67 
and declination (J2000) = —63? 73 with statistical un- 
certainty 0?06 (0?14) at the 68%(95%) confidence level, 
which is consistent with the position of PSR B1259— 63. 
The systematic uncertainty is estimated to be up to 
~40% (AbdoctaL 2009a). 

2.1. Light curve 

We then derived a 0.2-100 GeV light curve composed 
mostly of 3-day bins (Fig. [1]). Such bin size ensures 
enough photon statistics in each bin without loosing in- 
formation on short time-scale variability. Finer binning 
(e.g. 1-day) results in low significance for most of the 
days and high statistical uncertainties in deriving fluxes 
and photon indices, except for the flaring period that oc- 
cur between 2011 January 14 and February 3. For data 
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taken more than one month before and two months after 
the 2010 periastron passage, larger bins were employed 
for better visualization. From 2010 December 15 to 2011 
January 13, no evidence for 7-ra y emission was found, as 
first noted in IKong et al.l (j2011f) . We therefore derive an 
upper limit for the above quiescent period. Data points 
represent the flux values with TS>5, for which photon 
indices are also shown. Otherwise, an upper limit is de- 
rived. 

As shown in Fig. [TJ the 7-ray light curve is highly vari- 
able through the 2010 periastron passage: (1) The source 
started to be active in 7-rays about a month before 2010 
periastron passage (PI); (2) It remains undetected for 
about one month since mid-December (Ql); (3) Subse- 
quently, a major flaring period was identified; it peaks at 
~35 days after periastron. Having an average flux higher 
than that in PI by an order of magnitude, this flare lasted 
for only ~7 days (P2); (4) A second flare that peaks 
at ~46 days after periastron, however, lasted longer, so 
that the source was detected until end of February (P3 
and P4). To probe shorter time-scale variability during 
the flaring period (i.e., P2 and P3), we also produced 
the light curve with 1-day binning during January 14 to 
February 4 (Fig. [5]). Here 0.1-100 GeV photons were 
used to increase photon statistics. To better demon- 
strate the variability, we plot data points with TS>5 
(rather than the more standard criterium of TS>9) as 
flux values since a substantial number of data are of 
TS between 5 and 9. It can be seen that 7-rays from 
PSR B1259— 63/LS 2883 undergo rapid variations down 
to time scale of one day. Such variability is one of the 
fastest detected from any Galactic GeV source on the 
sky. No emission was detected during January 22-26, 
indicating that the source underwent a short quiescent 
period between the two major flares. 

We also show in Fig. [T]the orbital X-ray and TeV light 
curves. The GeV orbital light curve is different than the 
X-ray and TeV light curves, suggesting that the origin of 
GeV 7-rays is different than the others. 

2.2. Spectral analysis 

As the derived photon indices change significantly with 
time, we defined four periods that are indicated in Fig- 
ure [T] and Table [H We further divided the 0.1-300 GeV 
7-rays arriving during PI, P2, P3, and P4, respectively, 
into six energy bins of logarithmically equal bandwidths 
and reconstructed the flux using gtlike for each band in- 
dependently. A power-law (PL) model for each bin was 
assumed and the photon spectral index was fixed at a 
representative value F^ =2.8. While the emission is de- 
tected (i.e., TS>5) from 1.4 GeV to 20 GeV for PI, no 
7-ray source was apparent at the PSR B1259— 63 position 
in the four bins above 1.4 GcV (the derived TS values <5) 
in the likelihood analysis for P2, P3, and P4, indicating a 
cut-off at energy ~1 GeV during the flaring period. See 
Fig. 2] for the spectrum derived from a combined analy- 
sis of P2 and P3. We therefore attempted to fit the 0.1- 
100 GeV spectrum with a PL with an exponential cut-off 
(PLE), as well as a broken PL. We found that PLE de- 
scribes the spectrum even better during the periods P2, 
P3, and P4, by ATS=TSple-TSpl >8, i.e., >3cr in sig- 
nificance levels. The best-fit parameters are shown in Ta- 
ble [TJ In particular, the cut-off energies were found to be 
310±160 MeV, 550±330 MeV, and 250±95 MeV during 



3 



the periods P2, P3, and P4, respectively. On the other 
hand, the parameters provided by the broken PL with 
ah four parameters being free are not well constrained; 
we therefore do not consider t his model. 

At a distance 2.3 kpc (jNegueruela et al.l 1201 Ih . 
the average energy flux during the flares of 
~3xl0~^° erg cm~^ s^^ corresponds to the 7-ray lumi- 
nosity 1.9x10^^ erg s^^. Given the pulsar spin-down 
luminosity ^8x10'^^ erg s~^, the average 7-ray efficiency 
is about 25% during the flares. 

To demonstrate that GeV emission originates from a 
component that is temporally and spectrally different 
than the TeV emission, we put together the 100 MeV 
to 300 GeV spectra and the >300 GeV 7-ray spectra ob- 
tained for different orbital phases in Figure |4l It is clear 
that the GeV emission evolves differently compared to 
the TeV emission, assuming that the TeV behavior does 
not change dramatically between 2004 and 2010 perias- 
tron passages. 

2.3. Correlation between flux and photon index 

In Fig. [3] we plot 7-ray flux versus photon index, show- 
ing a correlation between these two quantities. We car- 
ried out a nonpar ametric correlation analysis. The com- 
puted Spearman rank correlation coefficient between two 
quantities is —0.7363. The probability that this coeffi- 
cient is different than zero is 0.9959. We have also calcu- 
lated the linear correlation coefficient (i.e., Pearson's r). 
This results in Pearson's r= —0.7874 and the probability 
that this coefficient is different than zero is 0.9986. 

3. DISCUSSION 

High-energy emissions from 7-ray binaries 
(PSR B1259-63, LS 5039, and LS I-h61°303) have 
been discussed using leptonic models (e.g., Tavani & 
Arons, 1997; Dubus 2006; Takata & Taam 2009) and 
hadronic models (e.g., Kawachi et al. 2004; Chernyakova 
et al. 2009). In leptonic models, PW particles (electrons 
and positrons) are accelerated at the shock where the dy- 
namical pressure of the PW and that of the stellar wind 
are in balance. These particles in turn emit non-thermal 
photons over a wide range of energies via synchrotron 
radiation (radio to GeV 7-rays) and inverse-Compton 
(IC) upscattering off star light (>10 GeV 7-rays), 
resulting in two peaks in the broadband spectrum. I n 
hadronic models (e.g., iChernvakova et aP I2006L 120091 ) . 
the inverse Compton emission of high-energy electrons 
resulting from Tr'^-decay is responsible for emission 
from optical up to TeV energies. In this case, the 
broadband spectrum is rather flat in the energy range 
of 0.1-100 GeV. 

During the pre-periastron epoch, GeV emission at a 
flux level ^ 6 X 10~^^ erg cm~^ s~^ is detected from 
mid-November through mid-December. The fitted spec- 
trum during this period (PI) is relatively hard (F^ ^ 2). 
Such hardness may be consistent with the predictions of 
both the hadronic and the leptonic models (see below 
for a discussion of the leptonic interpretation during this 
period). 

The onset of the flare-like GeV emission at the true 
anomaly 110° — 130° occurred close to the disk passage 
at the true anomaly 80° — 140°, suggesting that the origin 
of the 7-rays might be related to the disk passage. On 
one hand, the observed cut-off at several hundred MeV 



does not favor hadronic models. On the other hand, it 
may be difficult to explain the flare-like GeV emission us- 
ing simple leptonic models as well. First, leptonic mod- 
els predict a c ut-off around 100 MeV for a synchrotron 
spectrum fe.g.. iTakata fc: Taa m"2009^. Second, the stel- 
lar disk pushes the shock towards the pulsar when the 
pulsar encounters the disk. Although the increase of the 
magnetic field enhances the synchrotron power, it also 
reduces the synchrotron cooling time of the accelerated 
particles. As GeV photons are emitted by fast-cooling 
particles, these two effects should compensate each other. 
Consequently, the shift of the shock position due to the 
stellar disk cannot enhance the GeV emission to the level 
that we observe. 

Doppler boosting may provide a plausible mechanism 
to produce the GeV flares duriiig Jan uary 14 - Febru- 
ary 3, 2011 (see also lKong et al.l[2011[ ). Numerical sim- 
ulations in the hydrodynamic limit imply that the post 
shock bulk flow for the binary system can be acceler- 
ated into relativistic regime because of a rapid expansion 
of the flow in the downstream region (Bogovalov et al. 
2008). For example, Dubus et al. (2010) discussed the 
effects of Doppler boosting on the orbital modulations 
of the X-ray and TeV emission. If the bulk flow is rel- 
ativistic and oriented radially away from the star, the 
Doppler effect may be at work close to the true anomaly 
corresponding to the direction of the Earth 130°). 
This is approximately the time when the GeV flares 
were observed (i.e., 110° — 130°). The Doppler effect 
increases the photon energy as E = T>E' and the inten- 
sity as Ii, oc 2?'^+"/^, where 2? is the Doppler factor, and 
a — T — 1 the spectral index. Non-prime and prime nota- 
tions represent quantities in the observer and co-moving 
frames, respectively. 

In 7-rays, a-y ~ 1 — 2 is expected from synchrotron radi- 
ation models, which is consistent with the Fermi results. 
An enhancement factor of 5-10 in flux during flares (P2 
and P3) compared to the emission before periastron (PI; 
see table [T] for details) suggests V ~ 1.5-2. 

Suzaku obser vatio ns indi cate a low energy break, i.e., 
around 10 keV ^Uchivama et al.l l2009'). Because the pho- 
ton index of the synchrotron spectrum below the break is 
a = —1/3, the enhancement factor becomes D'^~^/'^ ~ 3. 
Thus, the enhancement in X-rays is suppressed compared 
to 7-rays. It is important to obtain simultaneous ob- 
servations in X-rays and gamma-rays during the flaring 
period to test the boosting model. 

In leptonic models, the correlation between F-y and flux 
(Fig. El may be understood as follows. If the Doppler 
effect is responsible for the flares, photons are boosted 
to higher energies and the observed synchrotron fiux is 
amplified. However, the Doppler boost does not affect 
much the photon fiux of the IC component, since more 
scatters between particles and incident photons would 
occur in the Klein-Nishina regime. Consequently, the in- 
crease in the IC flux by the Doppler boosting and the 
decrease due to the Klein-Nishina effect will compensate 
each other. The energy fiux of the synchrotron radia- 
tion becomes much higher than the IC radiation, caus- 
ing the 0.1-100 GeV spectrum to be dominated by the 
high-energy tail of the synchrotron spectrum. This gives 
rise to the steep spectrum (~3; corresponding to fiare 
emission) shown in Fig. [31 The emission before perias- 
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tron (with ~2) may be related to emission originat- 
ing just behind the shock where the bulk Lorentz factor 
is ~1. Some leptonic models (e.g. Takata et al. 2009; 
Dubus et al. 2010) have predicted similar energy fluxes of 
synchrotron and inverse- Compton radiation just behind 
the shock, hence r.y ~ 2, which is close to the observed 
value in the power-law fit. Although there is no evidence 
for two spectral components during the pre-periastron 
period (see Fig. |4]), this possibility cannot be excluded 
given the relatively low significance detection, i.e., ^5a, 
during this period. 

PSR B1259-63/LS 2883 is the third known 
binaries with significa nt detecti on in GeV, 
after LS I-F61°303 (lAbdo et all l2009bD and 
LS 5039 (|Abdo et al.l l2009c[ ). iT^ay be intuitive 
to compare the three 7-ray binaries: 

1. For LS 1+61 303 and PSR B1259-63, the GeV 
gamma-ray peak occurs after periastron (and be- 
fore apastron), though not in the same orbital 
phase. For LS 5039, the gamma-ray emission 
peaked at periastron; 

2. While no orbital phase- related spectral change has 
been reported in LS I+61°303, the "softer when 
brighter" behavior of GeV 7-rays has been found 
for both LS 5039 and PSR B1259-63; 

3. In all three systems, spectral cutoff is observed at 
GeV energies during at least part of the orbit; 

4. In all three systems, the X-ray and TeV 
light curves are more correlated whereas GeV 
7-rays always come out in different orbital 
phase s (|Anderhub et al.l 120091 : iTakahashi et al.l 
I20090 . This strongly suggests that GeV 7-rays orig- 
inate from components different than the other two 
energy bands. 

PSR B1259-63/LS 2883 is the only system among 
these three for which the nature of the compact object 
is certain. Given the similarities of some emission fea- 
tures of these systems, the results presented here should 
shed light on the GeV radiation mechanism of the other 
systems. 
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Table 1 

0.1-100 GcV 7-ray properties during different periods through the 2010 pcriastron passage. 



Period Year Date 


True Anomaly 


Mode 


1 Photon I'lux 




Energy Flux 


Photon Index 


Cutoll' Enerf 




ATS/ 












(erg cm ^ s 






(MeV) 


sig. 


sig. 


PI 2010 Nov 11 - Dec 14 


-115.3-2.1 


PL 


{9.9±4.1)xl0- 


B 


(5.8±1.6)xl0 


-11 


2.1±0.2 




31/5.60- 




P2 2011 Jan 14 - Jan 22 


111.7-120.0 


PL 


(1.5±0.2)xl0- 


B 


(2.9±0.3)xl0- 


-lU 


3.0±0.1 




163/12. So- 








PLE 


(1.3±0.2)xl0- 


6 


(2.7±0.3)xl0- 


-10 


1.8±0.6 


310±160 


184/13.6(7 


21/4.6(7 


P3 2011 Jan 26 - Feb 3 


122.2-127.2 


PL 


(1.4±0.2)xl0- 


B 


(3.1±0.4)xl0 


-lU 


2.8±0.1 




109/10.4(7 








PLE 


{1.3±0.2)xl0- 


6 


(2.8±0.3)xl0 


-10 


2.0±0.5 


550±330 


117/10.8(7 


8/2.8(7 


P4 2011 Feb 4 - Feb 21 


127.2-135.1 


PL 


(8.7±l.l)xl0- 


V 


(1.9±0.2)xl0 


-10 


2.8±0.1 




97/9.9(7 








PLE 


(7.5±l.l)xl0- 


7 


(1.7±0.2)xl0 
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Figure 1. X-ray, GeV and TeV observations of PSR B1259— 63/LS 2883. The error bars represent 68% confidence- level uncertainties, (a) 
Light curve derived from Fermi/LAT data taken from 2008 August 4 to 2011 February 21. For better visualization in the true anomaly 
scale, bins span variable time duration: between 2010 November 18 and 2011 February 9, three-day bins were used, apart from the quiescent 
period (Ql) between 2010 December 15 and 2011 January 13. Power law is assumed in deriving the flux with the photon index being 
free. For periods without significant detection (i.e., TS<5), 2 confidence-level upper limits were calculated assuming F^ = 2.1; (b) LAT 
photon index variation over time; (c) X-ray light curve over the 2010 and previous periastron passages. 2004 XMM-Newton observations 
(Chernyakova et al. 2006) are marked with triangles, 2007 Suzaku observations with circles, 2007 Swift observations with crosses, 2007 
Chandra observations with solid circles (Chernyakova 2009), and 2010/2011 Swift observations with stars (this work); (d) Evolution of 
the X-ray photon index. Notations arc the same as in panel (c); (e) Gamma-ray light curve above 1 TeV. The measurements were 
carried o ut by H.E.S.S. thr ough the 2004 (fluxes using a bin width of 2 days; Aharonian et al. 2005) and 2007 (monthly fiuxes) periastron 
passages ll Aharonian et aL|[2009i ). Also shown are upper limits derived from 2005 and 2006 observations IjAharonian et aLll2009l ). 
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Figure 2. Daily light curve from January 14 to February 4, covering flaring periods P2 and P3. For those days without significant 
detection (i.e., TS<5), two confidence-level upper limits were calculated assuming F-^ = 2.8. 
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Figure 3. Photon flux versus photon index in the 0.2—100 GeV band. Data are binned in the same manner as in Figure[Tl Data obtained 
during periods PI, P2, P3, and P4 are shown in black, red, blue, and green, respectively. 
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Figure 4. 100 MeV to ~3 TeV 7-ray spectra in dN/dE representation during different orbital phase periods. Data between 100 MeV 
and 300 GeV were obtained by Fermi/LAT through the 2010 periastron passage. Data above 300 GeV were obtained by the H.E.S.S. 
telescopes through the 2004 periastron passage (reproduced from Fig. 4 in Aharon ian et al. (2005)). The observations were carried out 
in the periods shown in each panel which are chosen to be in the similar orbital phases for the two energy bands in the same panel. The 
LAT data were divided into six energy bins of logarithmically equal bandwidths and the flux in each bin was reconstructed using gtlike 
independently, assuming power law within each bin and fixing F-y = 2.1 for PI and Ql, and F.^ = 2.9 for the flaring period (January 
14 - February 3). Two confidence- level upper limits were calculated for those bins with TS<5. Middle panel: 7-rays at the lowest bin 
100-380 MeV are (marginally) detected at TS=10, i.e. significance level ~3. Bottom panel: the lowest two energy bands are further divided 
into two bins each to increase the spectral resolution. 



